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ABSTRACT 

We present the first determination of the intrinsic three-dimensional shapes and the physical pa- 
rameters of both dark matter (DM) and intra-cluster medium (ICM) in a triaxial galaxy cluster. 
While most previous studies rely on the standard spherical modeling, our approach allows to infer 
the properties of the non-spherical intra-cluster gas distribution sitting in hydrostatic equilibrium 
within triaxial DM halos by combining X-ray, weak and strong lensing observations. We present an 
application of our method to the galaxy cluster MACS J1423. 8+2404. This source is an example of a 
well relaxed object with a unimodal mass distribution and we infer shape and physical properties of 
the ICM and the DM for this source. We found that this is a triaxial galaxy cluster with DM halo 
axial ratios 1.53 ± 0.15 and 1.44 ± 0.07 on the plane of the sky and along the line of sight, respec- 
tively. We show that accounting for the three-dimensional geometry allows to solve the long-standing 
discrepancy between galaxy cluster masses determined from X-ray and gravitational lensing observa- 
tions. We focus also on the determination of the inner slope of the DM density profile a, since the 
cuspiness of dark-matter density profiles in the central regions is one of the critical tests of the cold 
dark matter (CDM) paradigm for structure formation: we measure a = 0.94 ± 0.09 by accounting 
explicitly for the 3D structure for this cluster, a value which is close to the CDM predictions, while 
the standard spherical modeling leads to the biased value a = 1.24 ± 0.07. Our study proves that is 
not possible to disprove the manifestation of the DM with ~ 0.25 per cent of error of failing to reject 
the null hypothesis. Our findings provide further evidences that support the CDM scenario and open 
a new window in recovering the intrinsic shapes and desired physical parameters of galaxy clusters in 
a bias-free way. This has important consequences in using galaxy clusters as cosmological probes. 
Subject headings: galaxies: clusters: individual (MACS J1423. 8+2404), gravitational lensing, X-rays: 
galaxies: clusters 



1. INTRODUCTION 

Clusters of galaxies represent the largest virialized 
structures in the present universe, formed at relatively 
late times. They are an optimal place to test the pre- 
dictions of cosmological simulations regarding the mass 
profile of dark halos. In this respect their X-ray emission 
can be successfully used to constrain the mass profile in 
relaxed systems, where the emitt ing plasma is expected 
to be in hydrostatic equilibrium (|Sarazinlll988| ). 

Besides this method, gravitational lensing provides 
an opportunity to measure the cluster masses with- 
out invoking the assumption of hydrostatic equilib- 
rium implicit in the X-ray base d mass determinations 
(jMiralda-Escude fc Babull I1995I). Nevertheless, since 
lensing is sensitive to the integrated mass contrast along 
the line of sight, it is natural to expect mass overesti- 
mates due to fortuitous alignments with mass concen- 
trations which are not physically related to the galaxy 
cluster or departures of the DM halo from spherical sym- 
metry (|Gavazzill2005T ). 

So far most studies of both dark matter (DM) and 
intra-cluster medium (ICM) on galaxy clusters have been 
limited to the standard spherical geometry, i.e. the 
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length scales along the line of sight are taken as the spher- 
ical radii measured on the plane of the sky. First, the 
spatial resolution of Chandra and XMM-Newton X-ray 
satellites with their high sensitivity and large collecting 
area has allowed to resolve the core of the clusters, and 
have detected departures from isothermality and spheri- 
cal geometry of the ICM and DM: for example evidence 
for a flattened triaxial dark matter ha lo around five Abell 
galaxy clusters has been reported by iBuote fc Canizarei 
(1996). Second, numerical simulations predicts that DM 
halos show axis ratio s typically of the order of ~ 0.8 
l|Wang fc White) [20091 ) . disproving the spherical geome- 
try assumption. Third, since lensing is sensitive to the 
integrated mass contrast along the line of sight, depar- 
tures from the spherical assumption can justify the long- 
standing discrepancy between galaxy cluster masses de- 
termined fr om X-ray and strong gravitational lensing ob- 
servations ( Gavazzi I2005D . the latter being significantly 
higher than the former: since cluster mass measure- 
ments are sensitive to the assumptions about symmetry, 
this suggests that clusters with prominent strong lens- 
ing features are not spherically symmetric and prefer- 
entially elongated along the line of sight increasing the 
magnitude of the lensing. The galaxy cluster A1689 is 
a well-studie d example of a cluster show i ng this mass 
discrepancy (lAndersson fc M adeiski 2004; Lemz eet al.l 
2008; lRiemer-S0rensen et al.ll2009l: iPeng et al.ll2009D .~ 

While galaxy clusters are one of the key cosmological 
probes, this hinges on our ability to accurately determine 
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both mass and shapes of clusters. An accurate knowl- 
edge of the intrinsic cluster shape is required to constrain 
structure formation models via observations of clusters. 
Asphericity in the gas density distribution of clusters of 
galaxies is crucial in modeling X-ray morphologies and i n 
using clusters as cosmological tools (|Inagaki et al.lll995j ). 

Previous attempts to infer the elongation of the ICM in 
galaxy clusters through a comparison between X-ray and 
Sunyaev Zel'dovich (SZ ) observations (see, for example, 
Dc Filippis et al. 2005) relied on simple /3-model param- 
eterizations of the surface brightness profiles and on the 
isothermality of the same ICM. While these works sug- 
gested departures of the ICM from the standard spherical 
modeling, their prescriptions are subject to limitations 
implied in this modeling of the ICM. 

In the last few years the theoretical understanding 
of DM halos has also improved. While spherically av- 
eraged density profiles of DM ha los are well fitted b y 
a universal DM density profile (|Navarro et al.l 119971) . 
non-spherical effects on the mass function of dark halos 
have been studied in detail with numerical simulations 
(ILee fc Shanllarinlll^irTenkins et al.ll200ltl.Iing fc Sutol 
12002ft . In particular the work of iLee fc Sutol (|2003f ) de- 
rived analytical solutions for gas embedded in triaxial ha- 
los on basis of perturbation theory, and they recovered 
hydrostatic solutions for the gas density and tempera- 
ture profiles both for the isothermal and the polytropic 
equation of state. They also showed that the ICM and 
DM halos are well approximated by a sequence of con- 
centric triaxial distributions with different eccentricity 
ratio, which has been inferred analytically. 

Given that we observe just the two-dimensional pro- 
jected properties of galaxy clusters on the plane of the 
sky and any information along the line of sight is lost 
in the process of projection, it is in general impossible 
to derive the intrinsic three-dimensional shape of an as- 
tronomical object from a single observation. However, 
one can overcome the previous limitation by combin- 
ing observations in different wavelengths. In this pa- 
per we propose a new method where we compute gen- 
eral three-dimensional hydrostatic equilibrium solutions 
of the intra-cluster gas under the gravity of triaxial dark 
matter halos in a non-parametric way and by combining 
X-ray and lensing data. Indeed without any assumption 
for models on the gas density and (deprojected) tempera- 
ture profile, we assume that the triaxial DM density pro- 
file is well described by an analytical model (the general- 
ized Navarro, Frenk & White model, hereafter gNFW). 
Thanks to the results of numerical simulations, we know, 
indeed, sufficiently well the DM physics, which is in fact 
very simple, only depending on the gravity, unlike the 
physics of the baryons, which is also affected by sources 
of non-gravitational energy. Moreover we have removed 
the observational biases in the determination of the de- 
projected temperature (and consequently of the mass) 
by ad optin g the spectral-like temperature estimator (see 
Sect. I3.4[) . In this way we have a bias- free estimate of 
the deprojected temperature and, therefore, of the clus- 
ter mass. 

In the present paper we aim at describing both 
the theoretical aspects of the problem and astrophys- 
ical applications on a case study, the galaxy clus- 
ter MACS J1423.8+2404 (hereafter MACS J1423). This 
is a massive cool-core galaxy cluster at z = 0.539 



which does not show any significant sign of disturbance 
in its morphology and it appears to be very relaxed 
(lEbeling et al J 12(3011: iKartaltepe et all 120081 ). The study 
of lLimousin et al. (2009) proved that there is discrepancy 
between the 2D X-ray mass recovered under the assump- 
tion of the standard spherical geometry and the observed 
surface mass profile from lensing data, the former being 
systematically lower than the latter. In this perspective, 
motivated by the need to improve our knowledge of the 
3D physical properties of gas and DM in galaxy clusters 
and, especially, to understand how the 3D geometry will 
affect the estimated masses, we believe that MACS J1423 
is a optimal target for our analysis. In particular, we 
will combine X-ray, strong (SL) and weak (WL) lensing 
data in order to infer the geometry, physical parameters 
and distribution of the ICM and DM, which is a fun- 
damental question in bound systems (galaxies, galaxy 
clusters, dark matter halos) that form in an expanding 
universe, as well as to determine of the inner slope of the 
DM, since the cuspiness of dark-matter density profiles 
in the central regions is one of the critical tests of the 
cold dark matter (CDM) paradigm for structure forma- 
tion (|Navarro et al.lll997D . 

The paper is organized as follows. In Sect. [2] we 
present a model for the three-dimensional structure of 
galaxy clusters and how to constrain the 3D ICM and 
DM shape. In Sect. [3] we summarize the most relevant 
aspects of the X-ray data reduction procedure, we outline 
the method applied to determine the X-ray properties for 
MACS J1423, describing our spectral and spatial analy- 
sis. In Sect. 2]we apply a joint X-ray and lensing analysis 
in order to infer the physical properties of MACS J1423, 
while in Sect. [5] we discuss our findings. We leave to the 
appendix the discussion of some technical details of our 
data reduction procedure. 

Hereafter we have assumed a flat ACDM cosmology, 
with matter density parameter fiom = 0.3, cosmological 
constant density parameter £1a = 0.7, and Hubble con- 
stant Ho = 70km/s/Mpc. Unless otherwise stated, we 
estimated the errors at the 68.3 per cent confidence level. 

2. THREE-DIMENSIONAL STRUCTURE OF GALAXY 
CLUSTERS 

In order to infer the model parameters of both the ICM 
and of the underlying DM density profile, we perform a 
joint analysis for strong+weak lensing and X-ray data. 
Indeed the lensing and the X-ray emission both depends 
on the properties of the DM gravitational potential well, 
the former being a direct probe of the projected mass 
profile and the latter an indirect proxy of the mass pro- 
file through the hydrostatic equilibrium equation applied 
on the gas and temperature density. In this perspective 
we outline the methodology in order to infer physical 
properties in triaxial galaxy clusters. The general idea is 
straightforward: a) w e start with a triax ial DM density 
model as described in lJing fc Sutol ((2002) , which is rep- 
resentative of the total underlying mass distribution and 
depends on a few parameters to be determined, namely 
the concentration parameter c, the scale r adius r s and 
the inner sl ope of the DM a b) following ILee fc Sutol 
(2003l l200l . we recover the gravitational potential and 
surface mass profile k of a dark halo with such triaxial 
density profile c) we solve the hydrostatic equilibrium 
equation for the density of the ICM sitting in the gravi- 
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tational potential well previously calculated, in order to 
infer a theoretical three-dimensional temperature profile 
T gas in a non-parametric way 4) the joint comparison of 
T gas with the observed temperature and of k with the ob- 
served surface mass give us the parameters of the triaxial 
DM density model, and therefore all the desired physical 
properties of ICM and DM triaxial ellipsoids (see Fig. 
!)■ 

We start by describing in Sect. 12. li the adopted triaxial 
DM density and gravitational potential model, focusing 
on the relation between elongation of ICM and DM el- 
lipsoids. In Sect. 12.21 we present a toy problem in order 
to outline the effect of triaxiality on the physical observ- 
ables. 

2.1. Perturbative Expansion of the Triaxial Halo 
Potential 

In the present study, to parametrize the cluster mass 
distribution, we consider a tria xial generalized NFW 
model gNFW (j.Ting fc Sutdl2002h : 



p{R) = 



$cPc 



(R/R ) a (l + R/R ) 



3-a ' 



(1) 



where Rq is the scale length, S c is the dimension- 
less characteristic density contrast with respect to the 
critical density p c z of the universe at the present 
epoch, and a represents the inner slope of the den- 
sity profile; p CiZ = 3H(z) 2 /8itG is the critical den- 
sity of the universe at redshift z, H z = E z Hq, E z = 

[Q M (l + zf + (1 - fiji/ - 00(1 + z? + ^a] 1/2 , and 



Sr. = 



200 



3 ln(l + c) -c/(l + c) ' 



(2) 



where c = r 2 oo/r s is the concentration parameter, r s is 
the scale radius, x = r/r 2 oo- 

The radius R can be regarded as the the major axis 
length of the iso-density surfaces: 



R 2 



ir_ 

bl 



(a r > b r > c r ). (3) 



The ellipsoidal shape of a halo iso-density surface has 
been evaluated by defining the two eccentricities: 



e br = Wl- 



= \ 1- 



(4) 



and a r > b r > c r implies e& r < e Cr . The values of e 2 
(cr = b r ,c r ) measure the degree of the deviation of the 
ellipsoidal iso-density surfaces from the spherical ones 
along the corresponding principal axis direction. 

The gravitational potential of a dark halo with the tri- 
axial density profile ( eqn. [H can be written as a com- 
plex implicit integrals (|Binnev fc Tre maine 1987). While 
numerical integration is required in general to obtain 
the triaxial gravitational potential, small eccentricities 
(e 2 < e 2 <1) expected in cluster scale halos enable to 
approach this problem with the perturbative expansion 
and to infer analytical solutions. 

There are mainly two different ways to perform the el- 
lipsoidal perturbation, the equal-volume and the equal- 
length approach. In the former the perturbed ellip- 



soids have the same volumes as the unperturbed orig- 
inal spheres; in the latter the perturbed ellipsoids and 
the unperturbed spheres have the same length scales, i.e. 
the major (minor) axes of the perturbed ellipsoids will 
be coincident with the radius of the unperturbed sphere. 
iLee fc Sutol (|2003| ) developed the equal-length perturba- 
tion theory in order to provide an ellipsoidal correction 
to the standard spherical modeling of the observed clus- 
ters. While in the standard spherical model the length 
scale in the direction perpendicular to the line of sight is 
taken as the spherical tangential radii, i.e. measured on 
the plane of the sky, in their ellipsoidal model this tan- 
gential length scale is no longer the same as the line of 
sight length. Their equal-length perturbative expansion 
provides a framework which can be observationally more 
constrained than the equal-volume one because, for ex- 
ample, in SZ observations the cluster length scales in the 
direction perpendicular to the line of sight are measured 
(|De Filippis et al.|[200l . 

Her e we briefly summarize the findings of ILee fe Sutol 
( 2003) relevant for this study. Assuming a triaxial gNFW 
model for the DM (eqn. [TJ, they retrieved the following 
lst-order approximation for the gravitational potential 
$: 



$(u) wC Fxiu) + C e "^ eCr F 2 (u) 



+C 



e? sin 2 9 sin 2 d> + e 2 cos 2 I 



(5) 



-F 3 (u), 



where u = r/Ro, and C — 47rG(5 c p cr iti? 2 1 , and the three 
functions, F±(u), F 2 (u) 7 aTid Fs(u) are defined as 



F 1 (u) 



F 2 {u) 



a-2 



u .10 



2 1 



t + 1 



2-a 



dt 



(6) 



3 u Jo \t + 1 



2-a 



dt 



a-2 



F 3 (u) 



u 3 J 



dt. 



2-a 



(It 



(^+1) 



3-a 



(7) 



(8) 



In the previous equations CF\ (u) represents the spherical 
contribution to $(u), i.e., the gravitationa l potential for 
the case of the spher ical density profile (|Makino et al.1 
ll998tlSu"to et al.l["l998| ). CF 2 represents another spherical 
contribution due to the volume changes of the perturbed 
ellipsoidal density profiles from the spherical ones, while 
CF3 represents the non-spherical deviation of $(u) from 
the spherical potential and introduces a mild dependency 
of the potential (f> on 9. F3 being an order of magnitude 
smaller than F\ and F 2 , it can be neglected in most of 
t he applications. 

ILee fc Su to (2003) found that iso-potential surfaces of 
triaxial DM halos are still best well approximated as el- 
lipsoids and they model the iso-potential surfaces as tri- 
axial ellipsoids with the rescaled major axis length, £, 
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and the two eccentricities, tb T and e Cr : 



*(0 = c 



(9) 



and that the eccentricities of iso-potential surfaces e 2 are 
written in terms of their halo eccentricities e 2 as 



(a-2)F 3 (u) 



l-(a-2)F 1 (u) 



(10) 



Notice that e CT = eo-(e CT , u, a) unlike the constant for 
the adopted dark matter halo profile. In the whole range 
of u, e a /e a is less than unity (e a /e a ~ 0.7 at the center), 
and decreases mildly as u increases. This means that the 
intra-cluster gas is more spherical in the outer volumes 
than in the center, and it is altogether more spherical 
than the underlying DM halo. This is intuitively under- 
stood because the potential represents the overall average 
of the local density profile, and also because the gas pres- 
sure is isotropic unlike the possible anisotropic velocity 
ellipsoids for collisionless dark matter. 

Strictly speaking, equation ([TUf is valid only for e 2 <C 
1. This mathematical demand is just partially in agree- 
ment with the ellipticity expected for cluster scale ha- 
los: the intermediate-major and minor-intermediate axis 
ratios of halos are t ypically of the orde r of ~ 0.8 



(|Wang fc White! I2009D . In this perspective Eee\fcSuto| 
2003) provided the following fitting formula to infer the 



correct value e£ um of the eccentricity from the previous 
analytical expression e CT (eqn. [TO)) : 



1 + [0.1 + 0.05 log(l + it)]e^ 
[0.2 + 0.03 log(l + u)K r , 
1 + [0.1 + 0.09 log(l + u)]eg r 
[0.2 + 0.03 log(l + u)]e^ . 



(11) 



(12) 



Equations (TTT1) and ([12"]) are accurate within 10% errors 
for e a < 0.6, while within 20% errors for e a < 0.8. 

The iso-potential surfaces of the triaxial dark halo 
coincide also with the iso-density (pressure, tempera- 
ture) surfaces of the intra-cluster gas. This is sim- 
ply a direct consequence of the X-ray shape theorem 
(jBuote fc~C anizarcsl I1994T) : the hydrostatic equilibrium 
equation (fT5|) yields 



VP x V$ = Vp x V<f> = 0. 



(13) 



2.2. 



Constraining ICM and DM shape: a toy problem 

In order to understand how an elongation of the ICM 
and DM can be constrained by real observations, we 
present a simple toy problem and we analyze the effect 
of triaxiality on the physical observables. We consider 
an ideal ellipsoidal galaxy cluster with symmetry of rev- 
olution along the line of sight (i.e. = 0, e Cr + 1 0). 
The related gravitational potential can be obtained from 
equation 



*(u) « C 



Fi(u) ± e -^-{F 2 {u) +cos 2 9F 3 {u)} 



(14) 



where the positive and the negative signs in front of 
e 2 /2 correspond to the oblate [a r = b r > c r , e Cr = 

y/l — (c r /a r ) 2 )] and the prolate [a r < b r = c r , e Cr = 

yl — (or/br) 2 ] cases, respectively. Note that the z- 
direction is always chosen as the symmetric axis in the 
above expression. 

We consider an analytical model distribu- 
tion for the gas density in galaxy clusters, 
n e (r) = n (r/r c )~ e (l + r 2 /V 2 )~ 3/2/3+E/2 , with 
no, r c , e, /3 as free parameters. In order to repro- 
duce realistic galaxy clusters, we fix these parame- 
ters to the best fit values of the density profile of 
MACS J1423 under the assumption of spherical geome- 
try {n = 0.20cm- 3 ; r c = 21.7 kpc, e = 0.01, (3 = 0.55). 
We assume a spherical NFW model (i.e. a = 1 in eqn 
[T]) and fix the concentration parameter and the scale 
radius to c = 5.55 and r s = 253 kpc, which are the 
best expectation values we inferred for MACS J1423 in 
an X-ray-only analysis. We emphasize that we are not 
trying to reproduce exactly the physical properties of 
MACS J1423, but just a simple toy model for both the 
ICM and DM, whose parameters can be regarded as 
realistic for a typical massive galaxy cluster, in order 
understand how to constrain ICM and DM shape. 
We considered three different elongations of the DM 
halo (?7dm = {1,1-2,0.8}), corresponding to the case 
of spherical, prolate and oblate ellipsoid, respectively. 
From 77DM we inferred the corresponding elongation of 
the gravitational potential i]^ = ^(^dm, u, a), which 
is ranging from {1,1.09,0.90} to {1,1.03,0.96} moving 
from the center towards the X-ray boundary, and we re- 
calculated the corresponding gas density as n e oc i]^ ' 5 : 
indeed what is measured in observations is the emission 
measure K oc J n 2 dV, so it is pretty straightforward 
to recover the previous dependency of the gas density 
on the elongation along the line of sight (see also Ap- 
pendix A). In order to infer the dcprojcctcd temperature 
we solved the hydrostatic equilibrium equation (eqn 
[T5|) for the previous values of the elongation of both 
DM and ICM. Then we calculated the gravitational 
potential cp (eqn 1141) . the gas density n e , and the 
projected mass profile k for the same triaxial DM and 
ICM halos: the results of our analysis are reported in 
Fig. [TJ Concerning the temperature profiles compared 
to the spherical case we find an increase (decrease) AT 
of the temperature for prolate (oblate) ellipsoids, which 
is small but appreciable in the intermediate region 
(AT - 0.4 - 0.6 kcV for R ~ 200 - 600 kpc), but it is 
smaller near the center and the X-ray boundary. Bearing 
in mind that the shape of the temperature profile is the 
same as the underlying gravitational potential thanks 
to the X-ray shape theorem (eqn. I13p . the reasonable 
agreement among the temperature profiles for different 
values of t?dm in the inner and outer regions is intuitively 
understood because: a) the gravitational potential in 
the outskirts tends to be more spherical b) the tem- 
perature has been recovered by solving the hydrostatic 
equilibrium equation (eqn I15[) via integration of the 
derivative of the gravitational potential profile <p: it is 
possibly to show that for a large range of physical radii 
the relation T{R, r) BU )/T(R, t/dm = 1) =; ^(R = 0) 
roughly holds, therefore it is easy to understand why 
AT is larger in the central volumes than in the cooling 
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region, where the temperature drops of a factor of 2-3 
for a cool-core cluster. 




Fig. 1. — The radial profiles for the 3D temperature T gas (i?) in- 
ferred from the hydrostatic equilibrium equation (upper panel) and 
the projected mass profiles k (lower panel) for different elongation 
of the DM halo r?rjM = {!> 1-2,0.8} (points, dashed and dot-dashed 
line, respectively) . The distance of the points for the spherical case 
is representative of the true spatial resolution of X-ray and lensing 
data for MACS J1423. 



We point out that in Fig. Q]thc temperature and the 
projected mass profiles for the spherical case have been 
represented with points whose distance is representative 
of the true spatial resolution of X-ray and SL data for 
MACS J1423. Therefore we observe that the spatial res- 
olution is pretty good for both X-ray and lensing data, 
allowing to accurately measure effects of an elongation 
along the line of the sight. 

Concerning the projected mass profiles we observe that 
the dependency of the projected mass profile k on the 
ellipticity of the DM halo is quite straightforward, i.e. 
k(R) oc ?7dM) so the departures among different projected 
mass profiles as function of tidm are more pronounced in 
the central regions. 

Therefore from this toy model we find that the pro- 
jected mass profile near the center and the observed tem- 
perature profile in the intermediate regions are prone to 
the value of 77dm: so it easy to understand why both 
SL and X-ray data can jointly constrain the elongation 
of ICM and DM ellipsoids along the line of sight, being 
the former a proxy of the projected mass profiles in the 
internal regions and the latter of the temperature profile. 

As final test, we use this toy model in order to evalu- 
ate the systematics involved in our assumption that the 
triaxial ellipsoid is oriented (oblate or prolate) along the 
line of sight (see Sect. 13.4)) . We rotate the previous tree- 
dimcnsional ICM and DM ellipsoidal distribution of an 
angle of 30 degrees respect to the line of sight, we recover 
the desired two-dimensional physical properties and then 
we apply our approach. The systematics on the fitting 
parameters and on the axial ratios have been estimated 



<^ 3 per cent, which are within the statistical errors: in 
particular, the major-minor axial ratio and the scale ra- 
dius (the concentration parameter) get slightly underes- 
timated(overestimated). We conclude that this assump- 
tion does not strongly affect the estimated parameters 
a nd we point out that i t is ju stified in light of the results 
of lOguri fc Blandfordl (|2009f) . who showed that SL clus- 
ters with the largest Einstein radii constitute a highly 
biased population with major axes preferentially aligned 
with the line of sight increasing the magnitude of lensing. 

3. THE DATASET AND THE ANALYSIS 

We focus our analysis on the high-z cool-core galaxy 
cluster MACS J1423. 

Motivated by the need to improve our knowledge of 
the geometry, physical properties and distribution of 
the ICM and DM in galaxy clusters, we believe that 
MACSJ1423 is a optimal candidate for our study, be- 
cause it is very relaxed, with very good SL, WL and 
X-ray data. 

Here we briefly summarize the most relevant aspects 
of our data reduction and analysis. 

3.1. Strong and weak lensing analysis 

We refer to the findings of iLimousin et all (|2009l ) , who 
presented gravitational lensing and optical properties of 
MACS J1423 based on two bands of HST/ACS obser- 
vations complemented from ground by K band from 
the CFHT and Keck spectroscopy. They combined 
the strong lensing constraints from two spectroscopically 
confirmed multiply imaged systems with the weak lens- 
ing measurements from faint background galaxies via a 
parametric mass model. Their lensing analysis shows 
that the cluster has a unimodal mass distribution, sug- 
gesting that this an example of relaxed object. From the 
lensing analysis the cluster looks elongated with a minor- 
major axial ratio on the plane of the sky of 1.53 ± 0.15 
and position angle of 26 ± 2 degrees. The 2D projected 
mass map has been rebinned into elliptical annuli, whose 
ellipticity, ce ntroid and position ang le is the same as that 
inferred from lLimousin et al.1 ()2009f ) . Then we calculated 
average values of the elliptical symmetric projected mass 
profile k(R), R being the minor radius of the 2D ellipti- 
cal annuli, once we masked out the central 30 kpc, which 
is affected by the mass distribution of the cD galaxy. We 
also calculated the covariance matrix C' among all the 
measurements of k(R). 

3.2. X-ray data reduction 

We have analyzed a Chandra X-ray observation for 
MACS J1423 from the NASA HEASARC archive with 
a total exposure time of approx. 120 ks. We summa- 
rize here the most relevant aspects of the X-ray data 
reduction procedure for MACSJ1423. The observation 
(observation ID 4195) has been carried out using the 
Back Illuminated S3 chip of ACIS-S. We have repro- 
cessed the event 1 file retrieved from the Chandra archive 
with the CIAO software (version 4.1.2) distributed by 
the Chandra X-ray Observatory Center. We have run 
the tool aciss_proces_events to apply corrections for 
charge transfer inefficiency, re-computation of the events 
grade and flag background events associated with colli- 
sions on the detector of cosmic rays. We have considered 
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the gain file provided within CALDB (version 4.1.3) in 
this tool for the data in VFAINT mode. Then we have 
filtered the data to include the standard events grades 0, 
2, 3, 4 and 6 only, and therefore we have filtered for the 
Good Time Intervals (GTIs) supplied, which are con- 
tained in the f It 1. fits file. We checked for unusual 
background rates through the lc_sigma_clip, so we re- 
moved those points falling outside ±3cr from the mean 
value. Finally, wc filtered ACIS event files on energy se- 
lecting the range 300-9500 keV and on CCDs, so as to 
obtain an events 2 file. 

3.3. X-ray spatial and spectral analysis 

We outline the methodology of spatial and spectral 
analysis in triaxial galaxy clusters. The general idea is 
to measure the gas density profile in an non-parametric 
way from the surface brightness recovered by a spatial 
analysis, and to infer the observed projected temperature 
profile by a spectral analysis. 

The images have been extracted from the events 2 files 
in the energy range (0.5 — 5.0 keV), corrected by the 
exposure map to remove the vignetting effects, by mask- 
ing out the point sources. We determined the centroid 
(x c , y c ) of the surface brightness by locating the position 
where the X and Y derivatives go to zero, which is usu- 
ally a more robust determination than a center of mass or 
fitting a 2D Gaussian if the wings in one direction are af- 
fected by the presence of neighboring substructures. We 
constructed a set of n (n — 43) elliptical annuli of minor 
radius r m around the centroid of the surface brightness 
and with ellipticity €&< (r) fixed to that predicted from the 
cllip ticity e&'(r ) of the DM halo from SL data (see S ect. 
I3~4l Sect. Oand cqn. [TOj) . We discuss in Sect. 0T2] 
to which degree this assumption is consistent with the 
observed elongation of the X-ray isophotes. The minor 
radius of each annulus has been selected out to a maxi- 
mum distance i? S pat = 890 kpc, selecting the minor radii 
according to the following criteria: the number of net 
counts of photons from the source in the (0.5-5.0 keV) 
band is at least 200-1000 per annulus and the signal- 
to-noise ratio is always larger than 2. The background 
counts have been estimated from regions of the same ex- 
posure, which are free from source emissions. 

The spectral analysis has been performed by extract- 
ing the source spectra from n* {n* = 8) elliptical annuli 
of minor radius r* n around the centroid of the surface 
brightness and with ellipticity equal to that predicted 
from the ellipticity e' b (r) of the DM halo from SL data 
(see above). We have selected the minor radius of each 
annulus out to a maximum distance -R sp0 c = 880 kpc, 
according to the following criteria: the number of net 
counts of photons from the source in the band used for 
the spectral analysis is at least 2000 per annulus and cor- 
responds to a fraction of the total counts always larger 
than 30 per cent. 

The background spectra have been extracted from re- 
gions of the same exposure, and we have checked for sys- 
tematic errors due to possible source contamination of 
the background regions. This is done considering also the 
ACIS "blank-sky" background files: we have extracted 
the blank sky spectra from the same chip regions as the 
observed cluster spectra and scaled the blank sky spec- 
trum level to the corresponding observational spectrum 
in the 9-12 keV interval, where very little cluster emis- 



sion is expected. Then we have applied the aspect so- 
lution files of the observation to the background dataset 
by using reproject_events, so as to estimate the back- 
ground for our data. We have compared the two methods 
of background subtraction and we found that the differ- 
ence between the final results (e.g. the temperature) is 
negligible. In the following we have used only the results 
obtained using the local background. 

All the point sources have been masked out by both 
visual inspection and the tool celldetect, which pro- 
vide candidate point sources. Then we have calculated 
the redistribution matrix files (RMF) and the ancillary 
response files (ARF) for each annulus: in particular we 
have used the tools mkacisrmf to calculate the RMF, 
and the tool mkarf to derive the ARF of the regions. 

For each of the n* annuli the spec tra have been ana- 
lyzed by using the package XSPEC (I Arnaudl 119961 ver- 
sion 11.3.2) after grouping the photons into bins of 
20 counts per energy channel (using the task grppha 
from the FTOOLS software package) and applying the 
X 2 -statistics. The spectra are fitted with a single- 
temperature absor bed MEKAL model (Kaastra 1992; 
iLiedahl et al.l Il995h multiplied by a positi ve absorption 
edge as described in iVikhlinin et al.l (|2005|) . This proce- 
dure takes into account a correction to the effective area 
consisting in a 10 per cent decrement above 2.07 keV. 
The fit is performed in the energy range 0.6-7 keV (0.8-5 
keV for the outermost annulus only) by fixing the red- 
shift to the value obtained from optical spectroscopy and 
the absorbing equivalent hydrogen column density iVjj 
to the value of the Galac tic neutral hydrogen absor ption 
derived from radio data (jDickev fc Lockmanlll990l ). We 
consider three free parameters in the spectral analysis 
for the m— th annulus: the normalization of the thermal 
spectrum K m oc J dV, the emission-weighted temper- 
ature T* ro j m ; the metallicity Z m retrieved by employ- 
ing th e solar abundance ratios from IGrevesse fc Sauval 
(fTO) . 

The global (cooling-core corrected) temperature T ow 
has been derived in an elliptical region of ellipticity ey (r) 
and minor radius R, with 50 kpc < R < R S p C c, centered 
on the symmetrical center of the brightness distribution. 
We obtained a global cooling-core corrected temperature 
Tow = 6.34l°;^ keV and an abundance of 0.52 ± 0.07 
solar value for a reduced % 2 of 0.88 (356 degrees of free- 
dom). 

We cl assify this cluster as a st rong cooling core source 
(SCC) (jMorandi fc Ettoril[2007l ). i.e. the central cooling 
time icooi is significantly less than the age of the uni- 
verse t a ge,z at the cluster redshift (t coo i/t ag0 , z < 0.1): wc 
estimated a i coo i ~ 6 x 10 8 yr. As other SCC sources, 
MACSJ1423 show very low central temperature (~ 2 
keV) and a strong spike of luminosity in the brightness 
profile. The temperature profile is very regular, suggest- 
ing a relaxed dynamical state (see upper panel of Fig. 

02). 

3.4. X-ray spectral deprojection analysis in triaxial 
ellipsoids 

To measure the pressure and gravitating mass profiles 
in our clusters, we deproject the projected physical prop- 
erties obtained with the spectral and spatial analysis by 
using an updated and extended version of the technique 
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presented in iMorandi et al.l (|2007| ) . Here we summarize 
briefly the main characteristics of the adopted technique: 
(i) the electron density n e (r) is recovered both by depro- 
jecting the surface brightness profile and the spatially 
resolved spectral analysis obtaining a few tens of radial 
measurements in elliptical annuli; (ii) once a functional 
form of the DM density profile p — p(r, q) is assumed 
(i.e. a gNFW model, see eqn. [1]), where q = {q\, q%, ... qh) 
are free parameters of the DM analytical model, and the 
gas pressure Po at i? spC c is kept frozen to the observed 
value, the deprojected gas temperature T(q) is obtained 
by integration of the hydrostatic equilibrium equation: 



VP = -pV0 



(15) 



where p, = 0.6 is the average molecular weight, m p is the 
proton mass, P = pkT/pm p is the gas pressure, p is the 
gas density and (f> is the gravitational potential of the 
cluster. So T(q) = P(q)/n gas expressed in keV units. 

This approach is very powerful in order to infer the 3D 
temperature T(q), gas and DM mass profile, because it 
does not require any "real" spectral analysis for T(q), 
which could suffer of the poorness of the statistics and 
would need at least ~ 2000 net counts per annulus: we 
can determine the density in elliptical annuli even with 
very small counts (~ 200 — 1000). In other words we 
have an improvement (of about one order of magnitude) 
of the spatial resolution in the spectral analysis: indeed 
we remember that the total number n of elliptical annuli 
where we performed the spatial analysis is n = 43, while 
the total number n* of those where we performed the 
spectral analysis is n* —8. 

Given that in eqn. [3] the condition a r > b r > c r must 
always hold and the axial ratio of the elliptical DM halo 
are free parameters to be determined, i.e. unknown a 
priori and that might change orientation for each of the 
MCMC simulations (see below), we redefine new refer- 
ence axes (a',b',c'), which are simply azimuthally and 
longitudinally rotated of {0, ±90} degrees respect to the 
reference axes (a r , b r , c r ), and their orientation is frozen 
along the minor and major axis on the plane of the sky 
and the axis along the line of sight, respectively. Here- 
after, unless otherwise stated, we will always refer im- 
plicitly to the reference axes (a', b', c'), therefore ey (e&/) 
and e c i (e c / ) will refer to the eccentricity on the plane of 
the sky and along the line of sight, respectively. From 
Sect. I3.3l we remember that the ellipticity ey of the ellip- 
tical annuli where we performed the spectral and spatial 
analysis has been fixed to that predicted from the ellip- 
ticity eb' of the DM halo from SL data. Therefore, for 
the triaxial gNFW, we have: 



(c, r s ,a,e c >) 



(16) 



Concerning an X-ray-only analysis, the comparison of 
the observed projected temperature profile T* - 5 (Sect. 
3.31) with the deprojected T(q) inferred from hydrostatic 
equilibrium equation (eqn. I15[) . once the latter has been 
re-projected by correcting for the temp erature gradient 
along the line of sight as suggested in iMazzotta et al.l 
(2004), provides the probability distribution function of 
model parameters q via Markov Chain Monte Carlo 



e v = v/i- (f//«') 2 ; . 



•JX — (c'/a') 2 ) for oblate ellipsoids, 



(MCMC) algorithm, and therefore of T(q). The results 
for each of the MCMC simulations will be reported in the 
reference axes (a',6',c'). Even if an X-ray-only analysis 
might provide information on the desired model param- 
eters, given that we also behave lensing information, we 
will focus on a joint X-ray -(-lensing analysis (Sect. 14. ip . 

We point out that, in order to infer the electron density 
n e (r), we calculated the electron density n e (r) by depro- 
jecting the surface brightness profile in elliptical annuli 
of eccentricity £b'(r) and assuming first that an elonga- 
tion of the ICM ellipsoid halo along the line of sight is 
the same of the minor axis on the plane of the sky, and 
then we parametrize our ignorance about it through the 
following relation: 



n e (r) = h e (r) f] gaSlC '(r) 



(17) 



The proof of the previous eqn. [17] is reported in Ap- 
pendix A. 

4. APPLICATION OF OUR METHOD ON MACS J1423 
4.1. Joint X-ray+lensing analysis 

The lensing and the X-ray emission both depends on 
the properties of the DM gravitational potential well, the 
former being a direct probe of the projected mass pro- 
file and the latter an indirect proxy of the mass profile 
through the hydrostatic equilibrium equation applied on 
the gas temperature and density. In this sense, in order 
to infer the model parameters, we construct the likeli- 
hood performing a joint analysis for WL+SL and X-ray 
data, to constrain the properties of the model parame- 
ters q (eqn. [TBI of both the ICM and of the underlying 
DM density profile. 

We use the estimate of q from an X-ray-only analysis as 
proposal distribution to start a new MCMC simulation 
for the joint X-ray and lensing analysis: this makes the 
calculation faster, because we have a guessed starting 
value of q. 

The method works by constructing a joint X- 
ray+lensing likelihood: 



lens 



(18) 



being £ x and Acns the likelihoods coming from the X-ray 
and SL+WL data, respectively (see below). 
For £ x holds the following expression: 



a = 



exp{- X 2 /2} 
(27r)«*/2( CTl o 2 ...a n ,y 



(19) 



with x 2 equal to: 



E 

*=1 



CTproja(q) 



proj, 



(20) 



T* 



yT — (a'/c') 2 for prolate ones. 



being T* ro j j the observed projected temperature profile 
in the i — th ring and T 1 P roj,i( ( l) the convenient projection 
of theoretical 3D temperature Tj (q) in the j — th shell re- 
covered by applying the hydrostatic equilibrium equation 
(eqn. II 5|) and after correcting for the te mperature gradi- 
ent along the line of sight as suggested in IMazzotta et al.l 
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Fig. 2. — Example of a joint X-ray and lensing analysis (eqn. 
1 1 8 P> . In the upper panel we display the two quantities wh ich enter 
in the X-ray analysis spectral deprojection analysis (eqn. I20[l : the 
observed spectral projected temperature T* ro j m (big points with 

errorbars) and the theoretical projected temperature T pro j jIn (q) 
(diamonds). We also show the theoretical 3D temperature Tj(q) 
(points), which generates T pro i m (q) through convenient projection 
techniques. In the lower panel wc disp lay the two quantities which 
enter in the lensing analysis (eqn. 12 1 [ ): the observed surface mass 
profile fcj^ (points with errorbars) and the theoretical one fc m (q) 
(solid line). The distance of the points for both 7j(q) and is 
representative of the true spatial resolution of X-ray and lensing 
data, respectively. 



ter, for example the 3D temperature Tj = Tj(q) and 
the elongation e c /(e c ') of the ICM(DM) along the line 
of sight, just by relying on the hydrostatic equilibrium 
equation and on robust results of the hydrodynamical 
simulations of the DM profiles. In Fig. [2] we present 
an example of a joint analysis: notice that in the joint 
analysis both X-ray and lensing data are very well fit- 
ted by our model, with a total \ 2 = H>3 (15 degrees of 
freedom) . 

In fig [5] we present the joint probability distribution of 
c — r s and c — t)dm,c' , which represent marginal probabil- 
ity distributions of £(q, cr sys ). We observe that there is 
an anti-correlation between c and r s (j]DM,c')- Errors on 
the individual parameters (q, a sys ) have been evaluated 
by considering average value and absolute mean devia- 
tion on the marginal probability distributions of the same 
parameters. 





(|2004D . Aens reads: 

_ exp {-i[(fc 4 (q) - klfC-'Khiq) ~ K)]} 
Clcns ~ (27r) m */2|c*| 1 /2 ' 

(21) 

where C is the covariance matrix referred to the pro- 
jected mass profile from lensing data including system- 
atic effects (see below), | C| indicates the determinant of 
C, k* is the i — th observed measurement of the projected 
mass profile in the i — th elliptical annulus, fej(q) the the- 
oretical projected mass profile retrieved by our triaxial 
DM model (eqn. [1}, m* the total number of annuli. 

For the covariance matrix C, the below expression 
holds: 

C = C + C sys (22) 

being C the covariance matrix among the lensing mea- 
surements and C sys the covariance matrix arising from 
measurements of systematics, i.e. a bias parameter es- 
timator. We parametrize the systematics involved as 
C S ys = ofys I, where I is the identity matrix and <7 sy s a 
bias parameter estimator to be determined. er S ys repre- 
sents an example of latent variables (as opposed to ob- 
servable variables), i.e. variables that are not directly 
observed but are rather inferred (through a mathemat- 
ical model) from other variables that are observed and 
directly measured. We marginalized over (q, cr S ys) and 
therefore we have £ = £(q, a sys ). 

So we can determine the physical parameter of the clus- 
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Fig. 3. — Marginal probability distribution of c — r B and c — 
Vdm.c' (l e ft an d right panel, respectively). The solid(dashed) line 
represent the 1(2)-ct error region, while the big point represents the 
expectation value. 

Finally we computed the total mass enclosed in a ellip- 
soid £(Ra, eh', e c <) of major radius Ra and of ellipticitics 
e b >, e c i as M(q)(£(R A ,e b/ ,e c/ )) = J £ p toi (r, q) dV, where 
the radius Ra corresponds to a given overdensity A0: 
we considered the case where the overdensity is equal to 
2500. 

In table Q] we present the best model fit parameters 
for our analysis of MACS J1423. In Fig. H] we present 
an image of the core of MACS J1423 from optical (Hub- 
ble Space Telescope) observations, with overplotted the 
projected total mass contours computed from the gravi- 
tational lensing analysis (blue line) and from the X-ray 
surface brightness (green line). 

4.2. An independent test of our method: the elongation 
of the X-ray isophotes 

4 £(R&,e b i ,e c i) defines the ellipsoidal volume enclosing an av- 
erage density A times the critical density of the Universe. 
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TABLE 1 

Best model fit parameters of MACS J1423. The columns 1 — 5 refer to the best fit parameters c,r s , a,ri DMc > and a. 

WHILE THE LAST TWO COLUMNS REFER TO THE MASS AND RADIUS AT A = 2500, RESPECTIVELY. 



c r B a Vdm,c' CT sys M2500 -R2500 
(kpc) (gr/cm 2 ) (10 14 M Q ) (kpc) 



4.48 ±0.82 315 ±65 0.94 ±0.09 1.44 ± 0.07 0.006 ± 0.003 3.41 ±0.13 393 ± 12 



and the position angle of the major axis measured north 
through east in celestial coordinates can be estimated 
through the following equation: 



'ii = arctan ( — '- — ) + - (25) 



where A±(A + > A_) are the positive roots of the below 
quadratic equation: 

(M2,o-A 2 )( Aio ,2-A 2 )=/il 1 (26) 




Fig. 4.— Optical image of the core of MACS J1423 from F555W 
observation of the NASA/ESA Hubble Space Telescope, with over- 
plotted the projected total mass contours computed from the grav- 
itational lensing analysis (blue line) and from the X-ray surface 
brightness (green line). 
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In order to investigate the validity of our method we 
observe that the ICM eccentricity €(,' — £b'{&b' ,u,a) on 
the plane of the sky inferred from eqn. ([TO]) can be com- 
pared with the observed eccentricity of the X-ray sur- 
face brightn ess contours. We exploit th e iterate method 
proposed by lBuote fc Canizaresl (|1994l ) in order to mea- 
sure the flattening and orientation of the X-ray surface 
brightness. The parameters obtained from this method, 
rjM and 9m, computed within an elliptical region, pro- 
vide good estimates of the axial ratio T]m and the position 
angle 9m of an intrinsic elliptical distribution of constant 
shape and orientation. In order to determine these pa- 
rameters from an image of P pixels having rij counts in 
pixel i, one computes the moments: 

I p 

i=i 

where N = Y^h=i n ii an d ( x , V) is the centroid given by 
equation /ii — and /io.i = 0, respectively. Then the 
axial ratio reads: 

m - ^- (24) 



50 100 150 200 250 300 350 
R (kpc) 

Fig. 5. — The radial profiles for the axial ratio J? gas v on the plane 
of the sky predicted from eqn. HOI (solid line) and the observed 
elongation of the X-ray isophotes (points with errorbars). 



We compute rjM and 9m within elliptical annular aper- 
tures with increasing inner radii. Given that r\M is likely 
a function of the radial distance from the center, the el- 
liptical annuli should correspond more closely to a true 
isophote since only counts in the immediate vicinity of 
the isophote are used. We begin by defining a circular 
aperture for the first annulus {rjM = 1) about the cen- 
troid, which has been retrieved on a larger aperture, cov- 
ering the whole image. Then we compute the appropriate 
A*m,« for all the pixels in this elliptical aperture to obtain 
new values of em- Defining a new elliptical aperture with 
these parameters, we iterate until the parameters change 
by less than appropriate tolerances. These parameters 
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are also used to define the inner radius of the second 
elliptical annulus, for which we apply the same iterative 
procedure and so on: indeed we observe that the determi- 
nation of the desired parameters in each annulus affects 
the parameters in the more outer annuli. Characteriza- 
tion of the errors in this procedure has been performed 
via Monte Carlo randomization of uncertainties due to 
Poisson statistics of the surface brightness. We empha- 
size that the spread in value of tjm computed for many 
different t\m of the more inner annuli is a measure also 
of the importance of the systematic errors related to the 
chosen area of the annuli and departures of the bright- 
ness from the elliptical symmetry. Even if we did not 
find evidence for any position angle twists in the elliptical 
apertures, during the randomization of uncertainties that 
might occur because of the strong correlation among the 
parameters, so we decide to lock the position angle to the 
value inferred on a larger aperture which covers the whole 
image and randomize on this value (9m = 28.6 ± 1.2 
degrees), which is in good agreement with the value of 
iLimousin etafl (120091) (6& = 26 ± 2 deg r ees). 

As pointed out bv iBuote fc Canizaresl (|1994f ). for po- 



tentials whose shape changes with radius the bright- 
ness might not have exactly the same shape as the pro- 
jected potential, although they have the same three- 
dimensional shapes. Even if, strictly speaking, just for 
the case where <f> is stratified on concentric similar ellip- 
soids the brightness and the projection of cj> have exactly 
the same shapes, ind ependent of their th r ee-dim ensional 
radial distributions, IBuote fc Canizaresl (|1994f ) showed 
that also for small gradients in ellipticity, the projected 
shapes should closely approximate the similar ellipsoid 
case. This kind of systematics was implicitly taken into 
account in the previous estimate of the errors via Monte 
Carlo randomization. 

In Fig. [5] we report a comparison between the theoret- 
ical expectation for the axial ratio r7 gas ,b' on the plane of 
the sky from eqn. [10] and the observed elongation of the 
X-ray isophotes. The good agreement confirms the reli- 
ability of our method. We remember that our findings 
just rely on the hydrostatic equilibrium hypothesis and 
the assumption that the DM follows a triaxial gNFW, 
therefore the good agreement between observed and pre- 
dicted axial ratios indicates that the previous underly- 
ing assumption holds for this galaxy cluster. In other 
words, with our method, the desired physical properties 
are over-constrained by X-ray+lensing observations, pro- 
viding critical insights into our understanding of clusters, 
and critical tests of current models (i.e. DM halo model, 
hydrostatic equilibrium, 3D shapes) for galaxy clusters. 

We observe that in Fig. [5] a corresponding axial ratio 
7?DM,b' f° r the DM would be 1.53 ±0.15, value well above 
the axial ratio of the gas (?7gas,b' ~ 1-18 — 1.23). The 
probability that the DM axial ratio is in agreement with 
that of the gas has been discarded with ~ 0.25 per cent 
of error of failing to reject the null hypothesis via the x 2 - 
test, probing that the DM more elongated than the ICM 
(see Fig. |4]). The observed different elongation of halos 
of the visible matter (i.e. the ICM) constrained by X-ray 
data and of the gravitational one constrained by lensing 
observations adds strong evidence, in a very visual way, 
that the majority of the matter in the system is unseen 
and under the form of DM for general assumptions re- 



garding the behavior of gravity ( Cl owe et al.ll2006l) . 

5. DISCUSSION 

The main purposes of this work are to probe the shape 
of ICM and DM halos and value of the inner DM den- 
sity profile of MACS J1423, two quantities which are in- 
timately connected each other given their degeneracy, as 
well as to study the discrepancy between X-ray and lens- 
ing masses. Here we present the implications of our anal- 
ysis on the geometry, X-ray and lensing mass measure- 
ments (Sect. IBTTjl and the CDM scenario (Sect. IQ]) . 

5.1. Resolving the discrepancy between X-ray and 
lensing masses: probing the 3D geometry of ICM 
and DM 

Here we outline our findings in probing the 3D shape 
of ICM and DM, and the systematics involved in using a 
standard spherical modeling of galaxy clusters. In par- 
ticular, we will focus on the long-standing discrepancy 
between X-ray and lensing masses on clusters, showing 
that this is dispelled if we account explicitly for a triaxial 
geometry. 

MACS J1423 is a triaxial galaxy cluster with DM halo 
axial ratios r/DM,w = 1-53 ± 0.15 and ?7dm,c' = 1-44 ± 
0.07, where rjr>M,b' is the axial ratio of the DM on the 
plane of the sky inferred from lensing measurements, and 
?7dm,c' the axial ratio of the DM along the line of sight 
inferred through our analysis (see table [T]). Notice these 
elongations are statistically significant, i.e. it is possible 
to disprove the spherical geometry assumption. 

Here we are interested in emphasizing the importance 
of the 3D geometry on the desired physical parameters. 
If we assume a standard spherical modeling, we obtain 
a = 1.24 ± 0.07, c = 2.92 ± 0.77 and r s = 478 ± 136, 
values quit different from those in Table [1] emphasizing 
that the elongation of the source along the line of sight 
does affect the estimates of the physical properties: this 
proves that the systematics involved in neglecting elon- 
gation/flattening of the sources along the line of sight are 
important. The effects of geometry also may explain the 
large scatter of concentration para meter and inner sl ope 
of the DM found in the literature (|Oguri et al.ll2005H . 

We compare our concen tration parameter and sc ale ra- 
dius with the analysis of iSchmidt fc Alienl ((2007) , who 
used a circular NFW profile: they found r s — 170^20 
kpc and c = 8.27^1%. Again, this larger value of the 
concentration parameters is significant, showing that de- 
partures from sphericity are relevant in retrieving the 
desired physical parameters. In this perspective large 
concentration parameter values found in the literature 
for clusters with prominent strong lensing features can 
be explained by halos having the major ax is oriented to- 
ward the line of sight (jCorless et al.ll2009f) . 

Since lensing is sensitive to the integrated mass con- 
trast along the line of sight, departures from the spher- 
ical assumption could also justify the disagreement be- 
tween X -ray and lensin g mass profiles found in the lit- 
erature (Gavazzi 2005). Indeed, in Fig. [6] we compare 
the 2D mass enclosed within a circular apertures of ra- 
dius R for lensing, for an X-ray-only analysis under the 
assumption of spherical geometry, and from a joint X- 
ray+lensing analysis taking into account the 3D geom- 
etry. We emphasize the good agreement between the 
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masses inferred from lensing and a joint analysis based on 
triaxial modeling. On the contrary an X-ray-only analy- 
sis based on the standard spherical modeling clearly pre- 
dicts systematically lower masses by 20 -j- 35 per cent in 
the radial range 100 -j- 600 kpc. Notice that these dis- 
crepancies are increasing as we sample larger volumes of 
clusters: this is understood because of the larger value of 
r s (r s = 478 ± 136) for the spherical case respect to the 
triaxial one (r s = 315 ± 65). This confirms our insights 
about the role of the effects of geometry on the phys- 
ical properties and solve the long-standing discrepancy 
between X-ray and lensing masses of clusters. 

The application of our method on a larger sample of 
sources will allow to accurately retrieve cluster mass mea- 
surements, with important implications in the use of 
galaxy clusters as cosmological probes through the clus- 
ter mass function. As a note of caution, we remember 
that we focused on a relaxed cool-core object for which 
the underlying assumptions of our method, i.e. the hy- 
drostatic equilibrium hypothesis and that the DM follows 
a triaxial gNFW, are likely and has been proved to be 
tenable (Sect. I4.2|) . While we emphasize the importance 
of extending our method on other relaxed objects, nev- 
ertheless for non-relaxed galaxy clusters our approach 
might not improve significantly the knowledge of the 
physical parameters, being our underlying assumptions, 
in particular that the gas is in hydrostatic equilibrium, 
not wholly reliable for disturbed galaxy clusters. In this 
perspective, n-body simulations might address this issue 
for non-relaxed objects. 




200 400 

R (kpc) 

Fig. 6. — 2D masses enclosed within a circular aperture of radius 
R from lensing data (points with error bars), from an X-ray-only 
analysis under the assumption of spherical geometry (solid line 
with the 1-tr error hatched region), and from a joint X-ray+lensing 
analysis taking into account the 3D geometry (dashed line with the 
1-cr error gray shaded region). 



5.2. Probing the inner DM slope 

While it has long been recognized that N-body stud- 
ies of large-scale structure formation in cold dark matter 



(CDM) cosmologies form DM halos whose density pro- 
files are remar kably similar in shape over a wide range of 
halo mass (e.g. lDubinski fc Carlberdri99lHNavarro et al.l 
119971 ). a comprehensive physical explanation for the ori- 
gin of such a profile is still lacking; moreover, when it 
comes to the shape of the inner density profile of struc- 
t ures, no coherent pict ure has yet emerged. 

iNavarro et al.l (|1997l ) proved a general picture where 
this universal halo density distribution is character- 
ized by a relatively shallow powerlaw trend in the in- 
ner parts, pdm(t) oc r~ a , with a = 1. Subse- 
quently, seve ral studies suggested steeper central cusp 
(|Moore et al.l 119981 11999ft ; with one rece n t suit e of 
high-resolution simulations iDiemand et al.l (120041) ar - 
gue d for a = 1.16 ± 0.14 iTavlor fc Navarro! (|200l 
and IDehnen fc McLaughlin! (|2005|) developed a general 
framework where there is the possibility that a might 
be shallower than 1 at small radii, while other stud- 
ies suggested the lac k of cusp in the innermost regions 
(jNavarro et al.ll2004ft . 

On the observational side, efforts have been 
put on probing the central slope a of the un- 
derlying dar k matter distribution , through X-ray 
(lEttori et al.l 120021: lArabadiis et al.l 120021 ILewis et all 



20031: Zappacosta et al.ll2006ft : lensing ( Tyson et al 



Smith et al.l 120011: iDahle et all [20031: I Sand et al 



Gavazzi et all2003l: lGa"vazzi 2005; Sa nd et al. II 2004 12008 



1998; 



2002 



Brad adet al.l l2008t iLimousin et all 120081) or dynamics 



(jKelson et aillfooa iBrviano fc Saluccil 120061) . 



These 

studies lead to large scatter in the value of a from one 
cluster to another, but these determinations rely on the 
standard spherical modeling of galaxy clusters. Possi- 
ble elongation / flattening of the sources along the line of 
sight, as well as the degeneracy of a with other param- 
eters, i.e. c and r s , likely affect the estimated values of 
a. 

In this perspective, one of the main result of the pre- 
sented work is to measure a central slope of the DM 
a = 0.94 ±0.09 by accounting explicitly for the 3D struc- 
ture for MACS J1423: this valu e is c lose to the CDM 

1991 



predictions from INavarro et al.l (|1997ft (i.e. a = 1), but 
it is even in better agreement with the more recent nu- 
meral simulations of lMerritt et al.l i|2()()(i: , which predicts 
a slightly shallower inner slope. The value of the concen- 
tration parameter is 4.48 ± 0.82, in agreement with the 
th eoretical expectat ion from hydrodynamical simulations 
of lNeto et all (|2007l ) , where c ~ 4 at the redshift and for 
the virial mass of MACSJ1423, and with an intrinsic 
scatter of ~ 20 per cent. 

If we carry out a standard spherical modeling, we ob- 
tain the biased value a = 1.24 ± 0.07, value larger than 
that in table [TJ The different value of a in triaxial and 
spherical case shows that the systematics involved in ne- 
glecting elongation/flattening of the sources along the 
line of sight are relevant: this likely justifies the large 
scatter of a found in the literature. 

The lack of a flat core allows us to put a conservative 
upper limit on the dark matte r part icle scattering cross 
section: indeed lYoshida et all (|2000ft simulated cluster- 
sized halos and found that relatively small dark matter 
cross-sections (cdm = 0.1cm 2 g _1 ) are ruled out, pro- 
ducing a relatively large (40/i _1 kpc) cluster core, which 
is not observed in our case study. 
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We observe that the previous DM-only studies neglect 
the interplay between dark matter and baryons, which 
are present in observed structures. When it comes to 
account for the presence of baryons, the exact interplay 
between dark matter and baryon is not well understood, 
and different effects, i.e. cooling of the gas and dynamical 
friction acting on galaxies moving within the DM back- 
ground, compete in the final determination of the inner 
slope of the DM. While the former is expected to lead to 
a more concentrated dark matter density profiles via adi- 
abatic compression ([Gnedin et al.l l2004f) . the latter ca n 
heat up and soften the DM cusp (|E1-Zant et al.| [2004). 
In this analysis we did not subtract the mass contribu- 
tion from the X-ray gas to the total mass. Nevertheless 
the contribution of the gas to the total matter is small: 
the measured gas fraction is 0.06 — 0.07 in the spatial 
range 30 — 400 kpc, and the slope of the density profile 
is very similar to that of the DM beyond a characteris- 
tic scale r c ~ 20 — 30 kpc, a self-simil ar property of the 
gas c ommon to all the SCC sources (jMorandi fc Ettorj 
2007). This suggest that our assumption to model the 
total mass as a gNFW is reliable, and therefore a com- 
parison with DM-only studies is tenable. S i milar con- 
clusions have been reached by IBradac et al] ()2008l ) and 
iSommer-Larsen fc Limousin! ( 20091) . 

Both the observed different elongation of halos of the 
visible and gravitational matter (constrained by X-ray 
data and lensing observations, respectively) discussed in 
Sect. 14.21 and the value of the inner slope of the DM a 
strongly support the CDM scenario. 

6. SUMMARY AND CONCLUSIONS 

In this paper we have employed a triaxial halo model 
for the galaxy cluster MACS J1423 to extract more re- 
liable information on the three-dimensional shape and 
physical parameters, by combining X-ray and lensing 
measurements. We have obtained several significant re- 
sults. 

First, wc presented a new method in order to recon- 
struct the triaxial shape of both DM and ICM. It was 
demonstrated that the halo triaxiality can cause a signif- 
icant bias in estimating the desired physical parameters, 
i.e. concentration parameter c, inner slope of the DM 
a and total mass if a spherical halo model is a priori 
assumed for the model fitting. 

We focused on the implications of our method on the 
CDM scenario, proving that the value o f c and a are 
in ag reement with the CDM predictions (Navarro et al. 
once we properly accounted for the 3D shape of 



m agri 

Mm, 



the cluster. Departures of c and a from the theoreti- 
cal expectation of the CDM scenario found in the liter- 
ature can be explained by halos having the major axis 
preferentially oriented toward the line of sight. In par- 
ticular, accounting for the 3D geometry allows to resolve 
the long-standing discrepancy between X-ray and lensing 
masses in literature. 

Then we emphasized the implications of our analysis 
for DM. Our analysis proves that is not possible to dis- 
prove the manifestation of the DM with ~ 0.25 per cent 
of error of failing to reject the null hypothesis. 

This confirms our insights about the role of the ef- 
fects of geometry on the physical properties and allow 
to solve the long-standing discrepancy between X-ray 
and lensing masses of clusters. Nowadays the increas- 
ing precision of galaxy cluster observations makes the 
assumption of spherical geometry unlikely to be valid. 
The triaxiality of dark matter halos and ICM should 
systematically be taken into account for future analy- 
ses on galaxy clusters along the line we suggested. A 
relevant number of cosmological tests are today based 
on the knowledge of the mass and shape of galaxy 
clusters through X-ray measurements by assuming a 
spherical symmetry. Galaxy clusters play an impor- 
tant role in the determination o f cosmological par ame- 
ters such as the matter density (| Allen et al.l 120081 ). the 
amplitu de and slope of the density fl uctuations power 
spectr uiriVoevodkin fc V ikhlinin (|2004f ) , the Hubble con- 
stant (llnagaki et al.l Il995f) . to probe the nature of the 
dark energy (|Albrecht fc Bernstein! I2007t) and discrimi- 
nate between d ifferent cosmological sc enarios of struc- 
ture formation (jGastaldello et al.l 120071 ). It is therefore 
extremely important to provide a more general model- 
ing in order to properly determine the three-dimensional 
cluster shape and mass. The application of our method 
on a larger sample of sources will allow to accurately 
retrieve cluster mass and shape measurements, with im- 
portant implications for using clusters as cosmological 
tools. 
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APPENDIX 

APPENDIX MATERIAL: DEPROJECTING SPECTRA AND SURFACE BRIGHTNESS IN TRIAXIAL ELLIPSOIDS 

The deprojection technique decomposes the observed X-ray emission of the i-th annulus into the contributions from 
the volume fraction of the j-th spherical shells with j < i, by fixing the spectrum normalization (i.e. the emission 
measure Ki) of the outermost shell to the corresponding observed values. We can construct an upper triangular matrix 

V/, where the column vectors V 1 , V 2 , ...V™ represent the "effective" volumes, i.e. the volume of the j-th shell contained 
inside the i-th annulus (with j > i) and corrected by the gradient of n 2 , ^ inside the j-th shell fsee IMorandi et aT]|2007l 
for more details on the definition of effective volume instead of the more common geometrical one), so as: 

Ki oc n 2 c i dV = (v#n 2 ^j . (Al) 

In the previous equation n e = (n e .i, n e .2, n e ^ n ), being n* the total number of annuli, having internal (external) radius 
n n ,i , n ni 2 , •■• ,nn,n (?"out,i , T out ,2 , • r out,n); Ki is the XSPEC normalization of the spectrum in the i-th annulus; the 
operator # indicates the matrix product (rows by columns). Notice that the integral J\ >4 n 2 ^dV is of the order of the 
emission measure inside the i-th ring The inversion of this matrix allows us to determine n e ; . 

In the standard spherical modeling of the observed clusters the previous shells (annuli) and the matrix V\ refer to 
spherical annuli, while we are interested in extending the (de)-projection schema to concentric triaxial ellipsoids of 
reference axes a! , b' , c' with the axial ratios 77 gas ,b' = b' /a' and ?7gas,c' = c'/a' on the plane of the sky and along the line 
of sight, respectively. Instead of working on a orthogonal cartesian reference system (x, y, z) we chose another one, 
(x, y, z), such as x — x, y = j//?7gas.b'(f) and z = z/r] g&s , c >(r): in this new reference system it is intuitive to prove that 
the previous concentric triaxial ellipsoids can be regarded a spherical shells of radius a'(r). 

5 Hereafter we assume that the index j (i) indicates the shell (ring) of the source of radius (ri n ,r ou t). 
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Given the previous axis transformation, the relation between the effective volume matrices V and V in the reference 
system (x, y,z) and (x,y,z), respectively, reads: 

V = diag(ry gas , b 0#V#diag(?7 gas , c (A2) 

being diag(ry gaS!C /) a diagonal matrix with axial ratios (T? ga8 ,c',i) %as,c',2i % a s,c',n) on the main diagonal. Similar 
considerations hold for diag(?7g as ,b') 

Finally, in order to infer the gas density r? e in concentric triaxial ellipsoids of minor axis (a'i, a' 2 , a^), given the 
emission measure K in concentric elliptical annuli of minor axis (a^, a' 2 , a' n ), we have: 

K= (diag(T feaB)V )#V#diag(» 7gBBlC 0)#n= (A3) 

Notice that the matrix V refers to spherical annuli of radii (a[, a' 2 , a' n ). From eqn. IA3I it is easy to proof the 
dependency of the density on the elongation along the line of sight given by eqn. [TTj 
Similar considerations hold for the temperature (de)-projection. 



